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We investigated vortex dynamics in a pristine sample of type II superconduc- 
tor NbSe 2 by scanning tunneling microscopy at 4.2 K. The decay of the mag- 
netic field at ~ nT/s in our superconducting magnet induced the motion of 
vortices at ~ pm/s. Characteristics of the trajectories and speed variations of 
the moving vortices clearly revealed the potential landscape of vortices in the 
sample. Our vortex driving method may have wide application for the study 
of vortex matter on the nanometer scale by directly mapping an energy land- 
scape. 

Understanding vortex motion in superconductors plays a key role in developing useful 
superconductivity-based applications, and is crucial to answer many open questions about vor- 
tex matter. Although various techniques have yielded success in imaging vortices and have 
become instrumental in our microscopic understanding of vortex dynamics di]|2]|2]|4]l21), it is 
still challenging to capture the motion of individual vortices when inter- vortex spacing is less 
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than 100 nm and a few studies on this length scale have been reported to date. We used an 
extremely slow (~ nT/s) decay of a 0.5 T field (corresponding to inter-vortex spacing of about 
70 nm) in our superconducting magnet as the driving source for ultra-slow (~ pm/s) vortex 
motion in pristine NbSe 2 , and observed this motion by scanning tunneling microscopy (STM). 
We tracked vortices within a field of view, 400 x 400 nm 2 , over a week to extract highly time- 
and space-resolved data. Despite the use of pristine NbSe 2 , non-uniform motion was observed. 
The non-uniformity can be explained by the instability occurring when vortices exit the sample 
by overcoming surface barriers. This instability governed vortex motion allowed for mapping 
the potential landscape via trajectories and spatial variations of speeds of the moving vortices. 

When an external magnetic field between the lower and upper critical values is applied to a 
type II superconductor below its critical temperature, magnetic vortices self-arrange in the form 
of a triangular lattice, called a vortex lattice (VL), with its lattice constant © 



where B is the magnetic field penetrating the superconductor, and the magnetic flux quantum 
$ = 2.07 x 10~ 15 Tm 2 . We used the fi-dependance of the lattice constant of the VL to 
drive the motion of vortices. A small expansion of the lattice constant can lead to a large 
motion if the vortex system is comprised of many vortices. The magnetic field stored in our 
commercial superconducting magnet decays due to residual resistances at spot- welded joints in 
the superconducting coil. We found that the field decayed from 0.500 to 0.491 T over 21 days 
and 17 hours, which gives the average decay rate of AB/At ~ —4.2 nT/s —0.36 mT/day). 
Using a series inductor-resistor (L-R) circuit as a model (Fig.[T^), the magnetic field becomes 



where r = L/R. Eq. [2] with the initial/final magnetic fields measured and L = 12.4 H for our 
magnet, gives the decay time constant, r « 3.7 years, and the resistance, R — L/t ~ 0.1 (iVL. 
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This decay causes the triangular VL to expand (Fig. Q})), satisfying Eq. [TJ As a consequence 
of the continuity relation of vortex density and vortex flow, the speed of vortex motion may be 
calculated: 
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at a distance r away from the center of the sample. In this study, we used roughly a round 
disk of pristine NbSe 2 with a diameter of 5 mm (thickness = 0.5 mm). Since the STM tip was 
positioned near the center of the sample, a uniform motion, not exceeding ~ 10 pm/s, was 
expected. This model so far only predicts motion of a constant speed and straight trajectory. 
However, the motion could be affected by randomly distributed pinning sites (RDPS) (7) in the 
bulk or the surface barriers (SB) (8) along the edges. 

A home-built STM operating in vacuum at liquid helium temperature (4.2 K) was used. The 
sample was cleaved at room temperature under a pressure of ~ 6 x 1CT 7 mbar and transferred 
to the STM. When thermal equilibrium was reached, a magnetic field parallel to the c-axis of 
the NbSe 2 was raised from to 0.500 T at a rate of 0.100 T/min. The measurement was delayed 
4 hours after the initial field was stored in the superconducting magnet, to allow for transient 
relaxation of materials affecting the magnetic field. We operated the magnet in the persistent 
current mode during the measurement, subject to the slow decay. 

Vortices in an area of 400 x 400 nm 2 (128 x 128 pixels) were visualized in a dl/dV- 
spectroscopic image via a lock-in amplifier (Fig. [2k)- At a scan-speed of 600 nm/s (this gave a 
4 min completion time per frame), we did not observe significant distortions of the VL between 
adjacent frames, indicating that this scan- speed was sufficient to capture the motion of the 
VL. Over 17 days, we recorded consecutive spectroscopic images in the area, from which by 
automated data analysis we extracted the time-series of speed v(t) for all vortices ([9]). In this 
Letter, we present and discuss specifically the data taken during the last 7 days and ignore 
earlier data for two reasons. One is to avoid as much as possible any transient effects (|7Q|) 



after the initial ramp-up of the magnetic field. The other is that this was the longest continuous 
observation without any interruption of data collection, such as refilling the liquid helium dewar. 
This 7-day data series contains 2560 images of VL in total for a total number of v-data points 
of 102810. 

The gray-scaled image in Fig. [2k is the first spectroscopic image in the 7-day data series, 
showing a regular triangular VL with a lattice constant of 69 nm. As the magnetic field decayed, 
the VL moved mostly along one of its own principal axes (6Vl = 150°), while maintaining its 
ordered structure (see Supplementary Video 1 for the complete motion of VL over 7 days). As 
the time-colored trajectories of four vortices (labeled by 1, 2, 3, and 4) indicate, the motion 
of the VL was non-uniform in terms of speed and direction. Despite the non-uniform motion, 
the inter-vortex spacing remained almost constant throughout the measurement. Even at the 
highest velocity motion, 45 times faster than the averaged motion (marked as green circles in 
Fig. [2h and Figs. [3^,b), the ordered structure of the VL was still maintained. This rigidity of 
the VL indicate that a strong vortex-vortex interaction dominates distortions of the VL due to 
local pinning forces. After analyzing the 2560 images of the VL, we found that the lattice 
constant expanded by ~ 0.3 nm over 7 days, which is consistent with the field decay of ~ nT/s. 
However, we consider that this amount of change is negligible compared to the 69 nm lattice 
constant. Therefore, a single vortex motion represents the dynamics of the whole VL to a good 
approximation. This allowed us to combine all the time-series of individual vortices into a single 
time-series by sorting the observed speeds by time (black in Fig. [3k). Then we decomposed v 
into f par (red) and v per (blue) (v pai and f per are the speeds along and perpendicular to 6*vl 
respectively.). We found that (v), (u par ), and (v pei ) are 2.5, 0.9, and 0.1 pm/s respectively. 
The observed motion at ~ pm/s agrees well with the estimated speed based on the L-R model. 
However, to better understand the non-uniform motion of VL, it is necessary to include the 
effects of weak pining sites in bulk and surface barriers, as we will discuss later. 
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Strikingly, by overlaying the trajectories of all vortices over 7 days (color in Fig. [2b), we 
found that vortices repeatedly visited specific places (SPs) more often than other regions, al- 
though vortices arrived at the SPs along different paths. This characteristics also left the 2560- 
averaged VL image (gray- scaled), which represents the probability distribution of finding a vor- 
tex in space. The bright sites, where vortices were more likely to be found than in dark regions, 
clearly form a triangular lattice with its lattice constant close to that of the VL (Fig. [2^). 

We also investigated how the speed of the VL varied in space by defining the displacement 
of the VL, (X(t), Y(t)) (see Fig. [3b caption for the definitions). The (X(t), Y(t)) allows all 
trajectories to collapse into a single trace (discrete colors). By changing the variable from t to X 
in v pgx and v per (red and blue), one conveniently sees the variation of speeds in the (X, Y)-space. 
Down/up arrows mark the moments when vortices were at the SPs. When vortices were at the 
SPs, both v paT and v pcr were suppressed such that \J[v par (t)} 2 + [v per (t)] 2 < 3.5 pm/s. However, 
while vortices were transiting from SPs to SPs, accelerations/deccelerations occurred, showing 
fast motion such that \J[v pai (t)} 2 + [v per (t)] 2 > 5 pm/s. Therefore, the high probability and the 
suppressed speeds at the SPs indicate that the SPs are the minima of the potential landscape 
(PL) seen by the moving vortices, which has been preserved despite the continuous field decay 
over 7 days. 

In our observation of moving vortices in pristine NbSe 2 , there are two distinct features: 1) 
the non-uniform motion and 2) the manifestation/preservation of the PL. These two features 
could be caused by weak pinning sites, surface barriers, or the combination of both. If the 
weak pinning sites had played a dominant role in causing these two features, one would ob- 
serve collective flux creep motion (f5]), shortly after an external magnetic field is raised and set 
to a constant value. Even in a constant external magnetic field, the penetration of the mag- 
netic field into the superconductor is not instantaneous - the motion of vortices is impeded 
due to weak pinning sites, and is observable until they fully penetrate deep inside the super- 
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conductor. To check this, we ramped up the magnetic field from to 0.75 T, stabilized it by 
continuously supplying the current to the magnet in the persistent current mode (refer to the 
caption of Fig. for the detailed procedure), and immediately started the measurement, which 
continued over 2 days. One clearly sees quite reduced motion of vortices based on the overlaid 
trajectories (Fig.|4b). The (t> par ) and (f par ) were reduced by a factor of 10 compared to those of 
the moving VL. This suppressed vortex motion in a constant field right after the field ramp-up 
provides strong evidence that the two features were unlikely to be caused by weak pinning sites. 
Rather, this observation indicates that the crystal of NbSe2 we used in the study was a practi- 
cally pinning-free superconductor. Moreover, this experiment confirmed that the magnetic field 
decay was the main driving source of the vortex motion. 

If the influence of weak pinning sites cannot be significant, the most intriguing question 
is whether the surface barrier alone can cause the non-uniform motion and the manifesta- 
tion/preservation of the PL. Indeed the importance of the SB in vortex dynamics of high quality 
crystals of NbSe 2 has been emphasized in transport measurements (|771[72][73l) . Here we pro- 
pose a possible mechanism based on the interplay between the SB and the continuous magnetic 
field decay. As an external magnetic field decays, vortices exit the superconductor only when 
vortices overcome the SB along the edges (8). An exit event occurs at surface defects along the 
edges, as a bundle of vortices 5117751) . As soon as vortices exit, the vortices at the edges start 
re-arranging. This instability propagates through the superconductor via strong vortex-vortex 
interaction. As the field continues decaying, the instability keeps driving the vortex system with 
bundles of vortices leaving the sample at different locations based on local geometry. Therefore, 
the external force exerted on the vortex system chaotically varies in its magnitude, direction, and 
time. However, although a large number of vortices, 5N, were expected to leave the sample over 
7 days, it is still a small fraction (0.5 %) of the total number of vortices, N, when the VL image 
(Fig. [Hi) was taken; 5N = nr 2 sample 5B /$„ = 2.4 x 10 7 , and N = 7rr s 2 amplc 5 /$„ = 4.7 x 10 9 , 
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where r sample = 2.5 mm, 5B = 4.2 nT/s x 7 day « 2.5 mT, and 5 = 0.496 T. This negli- 
gible leakage of vortices should link to the preservation of the PL, as experimentally observed 
(Fig. |2})) . Therefore, each exit of a bundle of vortices can be considered as a perturbation to 
the vortex system. The continuing perturbations to the system cause the complicated patterns 
of trajectories and speed variation, but reveal the most energetically favorable configuration of 
vortices. 

For a full understanding of vortex dynamics, it is important to know the PL of vortex system. 
This is determined by many parameters - the superconducting properties (T c , H c ) of a super- 
conductor, the geometry, the details of pinning sites, temperature, and magnetic field. By this 
new method, the PL of pristine NbSe 2 in 0.5 T was clearly revealed for the first time by large 
statistical groups of vortices moving at very slow speeds. Our method will be useful in studying 
various vortex-pinning interactions - quenched disorders, periodic pinning sites, high-T c super- 
conductors - by directly mapping out an energy landscape, rather than by looking at instant 
motion of vortices. Moreover, in principle, the speed of the vortex motion can be controlled by 
adjusting the magnitude of the resistance R (Eq.[3]). This tunability of vortex motion combined 
with fast-scanning capability (1771 ITSl) will allow for direct systematic investigation of the dy- 
namic phases ([79]|20]) and phase transitions (|27][22]) of vortex matter in various superconductors 
on the nanometer scale. 
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Figure 1 : Model of driving source and vortex motion, a, L-R circuit corresponding to a field- 
decaying magnet due to non-zero R. L = 12 A H, R ps 0.1 iiVt, and r = L/R ~ 3.7 years (see 
text for the detail), b, Vortex motion model subject to the decay of a magnetic field. As a field 
B decays, vortices leave the sample. Observation within a field of view (rectangular area) is 
made at a distance r away from the center of the sample. The size of the field of view, vortices, 
and r are drawn for illustration, not to scale. 
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Figure 2: Trajectories of magnetically driven vortices in NbSe 2 . a, Collective and coherent 
motion of an ordered vortex lattice. The gray-scaled image, the first spectroscopic image of 
the 7-day observation, shows the configuration of the vortices in an area of 400 x 400 nm 2 
at 0.496 T. Subject to the magnetic field decay, all vortices moved together as the trajectories 
of four vortices (color) indicate. Vortices 1, 2, 3, and 4 remained within the field of view for 
5.1, 5.4, 4.5, and 4.1 days respectively. Trajectories are color-coded in time. The earliest time 
starts from red (1 st day), and ticks mark every 24 hour on the color bar. The 6>vl denotes the 
angle of one of principle axes of the VL. An ac modulation voltage AVbi as = 1 niV rms with 
/ = 1973 Hz, and a tunneling condition of J t = 0.1 nA and Vbi as = 3 mV between a STM tip 
and the sample were used, b, Overlaid trajectories (time-color scheme same to that of a) and 
the 2560-averaged spectroscopic image (gray-scaled) over 7 days. Notice that the places where 
the vortices more frequently passed by than other region form a triangular lattice. 
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Figure 3: Speeds of the VL motion, a, Speeds (v, t> par , and v pcr ) vs. time (t). t> par and v per 
are the components of the speed v of the VL, parallel and perpendicular to one of principal 
axes (6*vl = 150°) of the VL. The data of v, v par , and v pcr were smoothed by averaging over 
9 points to reduce noise, v and t> per are shifted vertically by ±10 pm/s respectively for clar- 
ity, b, Displacements and speeds. The displacements (X(t), Y(t)) of the VL are defined by 
X(t) = fo v pav (t)dt, and Y(t) = / * v pei (t)dt. The plot of (X(t), Y(t)) is discretely color-coded 
every 23 hour (= 1/ fo). The down/up arrows indicate when vortices were at the minimum of 
the potential landscape (bright sites in Fig. (2b). Dotted grid lines are drawn at every 69 nm, 
corresponding to the lattice constant of the VL at 0.496 T. The red and blue curves are the para- 
metric plots of (X(i), v pax (t)) and (X(i), v per (t)) respectively, where (X(t), v per (t)) is shifted 
vertically by —10 pm/s for clarity. 
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Figure 4: Turning off the vortex motion, a, Experimental procedure to maintain a con- 
stant magnetic field in a field-decaying superconducting magnet (Osheroff, D.; Private com- 
munication). First, while the heat switch is turned on (the residual resistance R becomes 
i?on ~ 20 fi), a current J flows through L, because the path across R on is now effectively 
open (r on = L/R on = 12.4 H/20 Vl ~ 0.5 s). Second, by turning off the heat switch to store 
the magnetic field B(I ) in the magnet, R on reduces to R g ~ 0.1 /xfi. In other words, r on 
increases to ~ 4 years (r D fj = 12.4 H/0.1 fxO,). Therefore, I from the external current source 
flows through R g instead of L. Lastly, one continues supplying I from the external current 
source to the magnet, even after the heat switch is off and B(Iq) is stored. Then two currents, 
one from the external current source and the other in the magnet, flow in the opposite direction 
with the same magnitude I across R Q g, which leads to quite a reduced power dissipation at 
R g, compared to that in the field-decaying case, b, Stationary vortices. According to the over- 
laid trajectories (color) over 2 days, vortices were prevented from moving while the procedure 
in a was employed. The continuous measurement time was significantly reduced, compared 
to the 7-day measurement. This is because the liquid helium boil-off rate increased with the 
magnet leads submerged in liquid helium to supply current during the measurement. The gray- 
scaled image was averaged over 704 VL images. The observation was made in an area of 
400 x 400 nm 2 under 0.75 T. The tunneling was achieved using I t — 0.1 nA and \4, ias = 3 mV, 
while scanning at a scan-speed of 551 nm/s. 
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